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Abstract We use the latest data to investigate observational constraints on the new 
generalized Chaplygin gas (NGCG) model. Using the Markov Chain Monte Carlo 
(MCMC) method, we constrain the NGCG model with the type la supernovae (SNe 
la) from Union2 set (557 data), the usual baryonic acoustic oscillation (BAO) obser- 
vation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 (DR7) 
galaxy sample, the cosmic microwave background (CMB) observation from the 7-year 
Wilkinson Microwave Anisotropy Probe (WMAP7) results, the newly revised H{z) 
data, as well as a value of 6*5 AO = 0.55) = (3.90 ± 0.38)° for the angular BAO 



rather stringent constraint. From the results, we can see a phantom model is slightly fa- 
vored and the probability that energy transfers from dark matter to dark energy is a little 
larger than the inverse. 

Key words: new generalized Chaplygin gas; angular BAO scale; cosmological observa- 
tions 

1 INTRODUCTION 

Cosmic observations suggest that the present universe is undergoing an accelerated state (Riess et al. 
1998, Perlmutter et al. 1999, Pope et al. 2004). In order to explain this, a component with negative 
pressure known as dark energy was proposed. The most simple and popular model for dark energy is 
cosmological constant (A). This model has successfully explained many phenomena while it indeed en- 
counters some theoretical problems. For example, the coincidence problem and the "fine-tuning" prob- 
lem. Therefore, many other models have been proposed including quintessence (Peebles et al. 1988a, 
Peebles et al. 1988b), holographic dark energy (Cohen et al. 1999, Li 2004), quintom (Guo et al. 2005, 
Feng et al. 2005), phantom (Caldwell 2002), brane world (DvaU et al. 2000, Zhu et al. 2005) and so 
on. Among these, the generalized Chaplygin gas (GCG) model as a unification of dark energy and dark 
matter is a good candidate (Bento et al. 2002a, Bento et al. 2002b). It has being widely studied and 
seems to be in agreement with different observational data (Zhu et al. 2004, Lu et al. 2009, Xu et al. 
2010a, Park et al. 2010). Moreover, GCG and the original Chaplygin gas (CG) model (Kamenshchik et 
al. 2001) can be connected with string and brane theory (Bilic et al. 2002). The equation of state (EOS) 
of GCG is expressed as 




A 
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Under Friedmann-Robertson-Walker (FRW) metric with energy momentum conservation equation , we 
can obtain (Amendola et al. 2003) 



Pgcg = Pgcg,o 



A + (1 - A^)a 



-3(1+q) 



(2) 



where A* = A/ Pq'qq q, pgcg,o is the energy density today. Besides the parameter ^, a is a new 
parameter in contrast to the original Chaplygin gas (CG) model. When a = 1, it becomes CG. The core 
of CG and GCG is when the scale factor is very small in the early universe, it acted as a dust, and it 
acts as the cosmological constant recently. The interacting CG model was discussed by Ref. (Zhang. H 
et al. 2006) and the interacting generalized Chaplygin gas model has been studied by Ref. (Setare et al. 
2007a, Setare et al. 2007b). Other models relevant to Chaplygin gas were proposed including modified 
Chaplygin gas (MCG) model (Chimento et al. 2004, Chimento et al. 2006), variable Chaplygin gas 
(VCG) model (Guo et al. 2007a), extended Chaplygin gas model (Meng et al. 2005) and new modified 
Chaplygin gas (NMCG) model (Chattopadhyay et al. 2008). 

Since the GCG can be equal to the interacting ACDM model (Fabris et al. 2004, Zhang. X et al. 
2006), a new generalized Chaplygin gas (NGCG) model which equals to a kind of interacting XCDM 
model was proposed as a unification of X-type dark energy with EOS parameter uix and cold dark 
matter (Zhang. X et al. 2006). The interacting new generalized Chaplygin gas model has been discussed 
by Ref. (Jamil et al. 2010). 

There are different kinds of observational data which can be used to constrain cosmological mod- 
els: type supernovae la, observational Hubble parameter data, baryon acoustic oscillation (BAO), cosmic 
microwave background data, lensing (Liao et al. 2012a) and so on. For BAO, using spectroscopic de- 
terminations of the redshifts of the galaxies, several detections at different redshifts have been studied 
(Percival et al. 2007, Percival et al. 2010). This method studies the three-dimensional averaged distance 
parameter Dy. The radial BAO Data was discussed by Ref. (Gaztanaga et al. 2009). Recently, a new 
determination of the BAO scale using the photometric sample of Luminous Red Galaxies in DR7 is per- 
formed (Crocce et al. 2011,Carneroetal. 2011). They get a value of 6'bao(^ 0.55) = (3.90±0.38)°, 
including systematic errors, for the angular BAO scale. It is the first direct measurement of the pure angu- 
lar BAO signal. Combined with previous BAO signals, it can break the degeneracies in the determination 
of model parameters. 

In this paper, we use the Markov chain Monte Carlo (MCMC) method to constrain NGCG model 
from the latest data including the BAO data at z — 0.55. In Section 2, we give a brief introduction of 
the NGCG model. In Section 3, we introduce the observational data. The full parameter space using 
different combinations of data is provided in Section 4. At last, we summarize our results in Section 5. 

2 THE NEW GENERALIZED CHAPLYGIN GAS MODEL 

We give a very brief introduction to the new generaUzed Chaplygin gas model in this section. For details 
of this model, please refer to Ref. (Zhang. X et al. 2006). Assuming the universe is flat with the FRW 
metric, the equation of state of NGCG fluid (Zhang. X et al. 2006) is 

A{a) 

Pngcg = — , (3) 

Pngcg 

where a is the scale factor and a is a parameter. NGCG fluid consists of dark energy px ^ a^'^'^^^™^-', 
where wx is the EOS parameter, and dark matter pdm ~ a '^- Naturally, the energy density of the 
NGCG can be considered as 



Pngcg 

with 



(4) 



A + B — Pngcg.o ' 
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where A and B positive constants. Considering the energy conservation equation, A{a) can be acquired 
as 

A{a) ^ -wxAa-^^^+"'''^^^+'''> . (6) 

We can easily see it becomes GCG when the EOS parameter of dark energy component wx is —1. On 
the other hand it becomes XCDM if a = 0. It is remarkable that a describes the interacting between 
dark energy and dark matter When a > 0, the energy transfers from dark matter to dark energy. On the 
contrary, the energy transfers from dark energy to dark matter in the case a < 0. Based on Ref. (Zhang. 
X et al. 2006), we take the radiation component into consideration, which dominated the early uni- 
verse. The whole density consists of NGCG component, the baryon matter component and the radiation 
component ptot = Pngcg + Pb + Pr- The Friedmann equation can be expressed as 



H{a) ^ HoEia), (7) 



where 



-3 



(1 -a"^'"^'') 



1 1/*? 



1 — fl/j — $1^ 

+nba-^ + 17,.a-^ . (8) 

The Ho, ilx, and Or represent Hubble constant, dimensionless dark energy, baryon matter and 
radiation today. Consist with the original paper, the parameter rj = 1 + a. 

3 CURRENT OBSERVATIONAL DATA 

Now, we introduce the methods of constraints on NGCG model by using the latest data. 

3.1 Baryon acoustic oscillation including data at z = 0.55 
For BAO, the distance scale is defined as (Eisenstein et al. 2005) 

and baryons were released from photons at the so-called drag epoch, the corresponding redshift Zd is 
give by 

where bi = 0.313(r2mo/i^)"° '^^^[l + O.607(r2mo/i^)°-^^1"^ and hi = 0.238(r2mo/i^)°-^^^ (Eisenstein 
et al. 1998). For observation of BAO, we choose the measurement of the distance radio (dz) at z = 0.2 
and z = 0.35 (Percival et al. 2010). The definition is given by 

dz= , s , (11) 
Dv{z) 

where rs{zd) is the comoving sound horizon. The SDSS data release 7 (DR7) galaxy sample gives the 
best-fit values of the data set (do. 2, do.35) (Percival 2010) 

do.2 \ ^( 0.1905 ± 0.0061 \ 

^0.35/ 0.1097 ± 0.0036 J ■ ^ ' 

The value of this BAO observation from SDSS DR7 can be calculated as (Percival et al. 2010) 

Xynatrix ~ matrix^ ma.tr ivL ^^matrix^ (13) 
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where the corresponding inverse covariance matrix is 



Cmatrix — ( i ■700-7 Q«n77 ) ■ (^^^ 



30124 -17227 \ 



-17227 86977 J 



Moreover, a new determination of BAO scale using photometric sample of Luminous Red Galaxies 
in the DR7 is presented (Crocce et al. 201 1, Carnero et al. 201 1). It make use of ^ 1.5 x 10^ luminous 
red galaxies with photometric redshifts. They get a value of Obao{z ~ 0.55) = (3.90 ± 0.38)° for 
the angular BAO scale including systematic errors. It is the first direct measurement of the pure angular 
BAO signal. The definition of Obao is expressed as (Sanchez et al. 2010) 

Obao = (15) 
where the comoving radial distance x(z) depending on NGCG model is defined as 

The corresponding of this BAO data is expressed as 

2 ^ [0BAO ~ 0.0681]^ 

Xangular 0.006632 ' ^ > 

where the unit has been transformed into radian. Since this data is independent of previous signals, the 
total can be acquired as 

XbAO Xmatrix ~^ Xangular' 

(18) 

3.2 Observational Hubble parameter data 

It is known that we use the distance scale of SNe la, CMB and BAO to constrain cosmological pa- 
rameters. However, the distance scale is determined by integrating the Hubble parameter, which can 
not reflect the fine structure of H{z). Thus, investigating the H{z) data directly reveals a more real 
scene. Many works have done using the newly revised H{z) (Xu et al. 2010b, Ma et al. 2011). The 
measurement of Hubble parameter data as a function of redshift z depends on the differential ages of 
red-envelope galaxies 

\ dz 
1 + z at 

We choose 12 data from this method (Riess et al. 2009, Stern et al. 2009). In addition, the data can be 
obtained from the BAO scale as a standard ruler in the radial direction. So, we choose another three data 
at different redshifts (Gaztanaga et al. 2008). 

The value of the H{z) data can be expressed as 

,2 _ [H{z,)-Hobs{z,)]^ 

Xh — 2^ 2 ' ^'^^) 

i=l '^i 

where is the la uncertainty of the observational H{z) data. 

3.3 Type la supernovae 

SNe la has been playing an important role in studying cosmology since it first revealed the accelerated 
expansion of the universe. The recent data (Union2) is given by the Supernova Cosmology Project (SCP) 
collaboration including 557 samples (Amanullah et al. 2010). The new data has been used to constrain 
cosmological models (Xu et al. 2010c, Xu et al. 2010d, Liao et al. 2012b). In practice, people use 
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the distance modules of supernovae to reflect the cosmological model and constrain the cosmological 
parameters. The distance modules is determined by the luminosity distance 

M = 51og(dL/Mpc) + 25, (21) 

where is the luminosity distance. In the flat universe, it is connected with redshift which is a obser- 
vational quantity 

dL = {l + z) [ dz'/H[z'). (22) 

JO 

We choose the marginalized nuisance parameter (Nesseris et al. 2005) for ~x^: 

Xm.^A-—, (23) 

where A = ^ (^'^^''^ - ^''^""y) Vcrf^ B = Y.f' (/x^^'^ - /^*'^<=°'y)/cr2, Q = ^f" 1/crf, a, is the 
la uncertainty of the observational data. 

3.4 Cosmic microwave background 

For CMB, the acoustic scale is related to the distance ratio and is expressed as 

la = TT '-^ , (24) 

where rs{z^) = Ho~^ Cs{z)/E{z)dz is the comoving sound horizon at photo-decoupling epoch. 
The CMB shift parameter R is expressed as (Bond et al. 1997) 



dz 



(25) 



/o Eiz). 

where the redshift z^ corresponding to the decoupling epoch of photons is given by (Hu et al. 1996) 

z, = 1048[1 + 0.00124(f]b/i^)-°-^3«(l + giin^^oh^^)], (26) 

where .gi = 0.0783(r!6/i2)-o.238(i ^ 39.5(r!fc;i2)-o.763)-i^ = 0.560(1 + 21.1(r!fc/i2)i-8i)-i. 

For the CMB observation, we choose the data set including the the acoustic scale (la), the shift 
parameter (R), and the redshift of recombination (z*). The WMAP7 measurement gives the best-fit 
values of the data set (Komatsu et al. 2011) 

fla\ ( 302.09 ± 0.76 \ 
PcMB = M = 1-725 ±0.018 . (27) 
\z,J \ 1091.3 ±0.91 / 

The value of the CMB observation can be calculated as (Komatsu et al. 201 1) 

XcMB = ^PcmbCcmb ^APcmb, (28) 
where APcmb = Pcmb ~ Pcmb, and the corresponding inverse covariance matrix is 

2.305 29.698 -1.333 
Ccmb"^ = I 29.698 6825.270 -113.180 | . (29) 
-1.333 -113.180 3.414 
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The NGCG Model 



SNc+CMB+BAO+H ( z ) 



"■'^^■^*_0.0010 ■I -0.0014 
"•J-J-'^y_0.0075^-"'J-0.0102^'''^-' 

^■"-^ ' -0. 0502 -0. 0716 

0.7297i---(l.)l«;»-(2.) 

13.73l«-(l.)l-?°(2.) 

o.2703;S;»=^^,(i.)+S;°«J{2.) 

71.03l^3f,(l.)l^,^^(2a) 



SNc+BAO+CMB 
0.02241-;- (1.)+--- (2.) 

o-"**^roS8?(i'^)lS:Sif8(2-) 

-1 0297+''-l''''^ri<Tl+°-^*^*r2tTl 
1 0076+''-''^''2ao-)+''-''''31('2cT'l 
O.T232+°;«3-(l.)+°---(2.) 

i3.75+S;f,{i.)+S;-(2.) 

0.2768+"-°2^»{la)+0-'"*« 



CMB+BAO+H(z) 

0225+" ™^^ flCT'l+'J-'"'^'*(2CT) 
"■"^^''-0.0011 '■-'''•'-O.OOIS^^'''' 



fib ft 

1 
fix 

Age/Gyr 

n„, 

ffo 



70.29 



-0.0305 ^ 
+ 3.86 



(la) 



-0.0408 
+ 6.61 



(2cr) 



(2<t) 



-0.0097'- ^-0.0138 
1 11t;fi+0-2786 ,, %+0.3588 /,„> 
-1.1156_Q 3jj,^{l<T)_p^^g5{2<T) 

1 0170+'''°*''''(lal+°-°'"'3(2<T'l 
-^•"-^ "J_0. 0598 >■ •'-0.0907 

0-7367_g Q3g2(l<T)_^ „5^3{2<t) 
13.70;0'^l(l<T)+03|{2a) 

0.2633;°;Sf,?(i.)+S;S=«{2.) 

72.36t^ l«(l.)r 



Table 1 The best-fit values of parameters r2b^^, f^DM^^ ^^x, ?7, f^x, Age/Gyr, ri„i, and 
for NGCG model with the l-cr and 2-(t uncertainties, as well as Xmin' '^he data sets 
SNe+CMB+BAO+iJ(z), SNe+BAO+CMB, and CMB+BAO+i7(z), respectively. 



4 CONSTRAINTS ON NGCG MODEL VIA MCMC METHOD 

In order to constrain the parameters of NGCG model. We take the usual maximum Ukehhood method 
of fitting with the Markov Chain Monte Carlo (MCMC) method. We take the Metropolis-Hastings 
algorithm with uniform prior probability distribution. By using Monto Carlo method, we generate a 
chain of sample points distributed in the parameter space according to the posterior probability, then 
repeat the process until the established convergence accuracy is acquired. In our testing, the convergence 
of the chains i? — 1 is set to be less than 0.003 which is small enough to guarantee the accuracy. The code 
we use is based on CosmoMCMC (Lewis et al. 2002). We combine the SNe, BAO, CMB and B{z) data 
by multiplying the likelihood functions to constrain the NGCG model. The total -y^ can be expressed as 



XsNo 



Xbao + XCMB + xlf- 



(30) 



We show the 1-D probability of each parameter (fJb^^, <-^x, ^x, Age/Gyr, fim, 

Hq) (fib, ^dm, and ilm are responding to baryon matter, dark matter, dark energy and the 
total matter) and 2-D plots for parameters between each other for the NGCG model with SNe 
+ CMB + BAO + H{z) in Fig. [T] The best-fit values of NGCG model parameters with the 
four kinds of data above are cJx=-L0510iS; }i^(lcr)lg-^^^^(2cr), r;=L0117ij] E]^Jj^(lcr)i;] J]??^(2cr), 

f],/i2=0.0224tO:0°}^(la)+0:^o}I(2a), 17.^=0.7297tror7^(l-)t°o;Sl^l(2^), and the xLi 543.668. 
We also constrain NGCG model with other combinations of data for comparison: SNe + CMB + BAO 
and CMB + BAO + H{z). They are showed in Fig.|2]and Fig. [3] The best-fit values of each parameter 
with the l-cr and 2-a uncertainties and the Xmin are presented in Table 1. We can see the combination of 
the latest data gives a rather tight constraint on NGCG parameters especially for 77 which describes the 
interacting between dark energy and dark matter The best-fit of ujx is slightly smaller than -1, which 
makes a little difference from Zhang et al. 2006, and ry is a little larger than 1 . Considering the uncer- 
tainty, these indicate the scene of the universe that the NGCG model reflects is near the ACDM model. 
On the other hand, the results still accommodate a interacting XCDM model. In 1-a range for SNe + 
CMB H- BAO H- = (-1.2195, -0.8947), a = - 1 = (-0.0385, 0.0585). The resuh of 

indicates the dark energy acts as a phantom is slightly favored. The constraint on a indicates the proba- 
bility that energy transfers from dark matter to dark energy is a little larger than the inverse. Moreover, 
since a is constrained to a very small range near 0, the interacting between dark sectors seems very 
weak. The results also totally rule out the original CG model and the VCG model which require a = 1. 
This is in agreement with (Zhang. X et al. 2006, Wu et al. 2007). 
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Fig. 1 The 2-D regions and 1-D marginalized distribution with the 1-cr and 2-cr contours of 
parameters n^h^, JIdm^i^ ^x^ V' ^x, Age/Gyr, fi^j and Hq in NGCG model, for the data 
sets SNe+CMB+BAO+H(z). 



5 CONCLUSION 

With the MCMC method, we give constraints on the NGCG model from the latest data including the 
SNe la data from Union2, BAO data from SDSS DR7 galaxy sample, and CMB observation from 
WMAP7, as well as the newly obtained angular BAO signal OBAoiz = 0.55) = (3.90 ± 0.38)°, 
which is the first direct measurement of the pure angular BAO signal. The best-fit values of the 
parameters for NGCG are wx=-1.0510l°-i^^^(lcr)t°-^226(2cr), 77=1.0117t°-°^^^(la)tEJ°??^(2(T), 

17,;^2=0.0224lSJ:™12o(la)t°:E5;;}I(2a), Ox=0:7297iro2f6(l'^)tS;S^o^(2^)- and the x?„i„ =' 543.668, 
which give a quite tight constraint. From the results we can see ACDM model is near the best-fit point, 
it remains a good choice for explaining the observation. However, the NGCG model permits a interact- 
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Fig. 2 The 2-D regions and 1-D marginalized distribution with the 1-cr and 2-cr contours of 
parameters n^h^, JIdm^^ wjc, rj, ilx, Age/Gyr, fi^j and Hq in NGCG model, for the data 
sets SNe+CMB+BAO. 



ing XCDM model. The constraint results seem to slightly favor a phantom and the energy transfers from 
dark matter to dark energy. Moreover, since a is constrained to a very small range near 0, the interacting 
between dark sectors seems very weak. The results are consist with the situations in Ref. (Guo 2007b et 
al. , Chen 2010 et al. , Feng 2008 et al. , Cao 2011a et al. , Cao 201 lb et al. ). These papers as well as 
ours manifest current observations can not distinguish the directions of energy transforming well, thus 
the coincidence problem can not be solved so far in this way. For future study on coincidence problem, 
we hope more data and other independent observations can give the distinguishability. Besides, since 
the constraint results for these interacting XCDM models slightly favor the energy transfers from dark 
matter to dark energy, we may suspect whether this way is valid for explaining the coincidence problem. 
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Fig. 3 The 2-D regions and 1-D marginalized distribution with the l-cr and 2-a contours of 
parameters Vl^h?, SIdm^^ wx, r], Q,x, Age/Gyr, 0,^, and Hq in NGCG model, for the data 
sets CMB+BAO+ff (z). 
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